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Abstract. - We present a self-consistent theory of Anderson localization that yields a simple 
algorithm to obtain typical local density of states as an order parameter, thereby reproducing 
the essential features of a phase-diagram of localization-delocalization quantum phase transition 
in the standard lattice models of disordered electron problem. Due to the local character of our 
theory, it can easily be combined with dynamical mean-field approaches to strongly correlated 
electrons, thus opening an attractive avenue for a genuine non-perturbative treatment of the 
interplay of strong interactions and strong disorder. 



After more than four decades of vigorous efforts, there remains little doubt that disorder- 
driven metal-insulator transitions (MITs) bear many similarities to more familiar critical 
phenomena. The basic physical process involved was identified by Anderson pp, who first 
discussed the localization of electronic wave functions as a driving force behind such MITs. 
Further theoretical progress has been slow, partly due to ambiguities in identifying an appro- 
priate order parameter for Anderson localization. Nevertheless, important information was 
obtained by using scaling approaches based on 2 + e expansions [2] , which were subsequently 
extended [3] to incorporate the interaction effects. 

There are several reasons why the existing theories remain unsatisfactory. Most impor- 
tantly, the MITs generically take place at strong disorder where the energy scales associated 
with both disorder and the interactions are comparable to the Fermi energy, in contrast to 
what happens in perturbative 2 + e-expansion approaches. As a result, well defined precursors 
of the MITs are seen even at very high temperatures, as experimentally demonstrated in many 
systems These features include not only the scaling behavior of various quantities, but 

also the breakdown of the Matthiessen rule and the Mooij correlation 0|. To understand 
such global behavior, a mean-field like formulation would be advantageous, but it should be 
one that can incorporate both the strong disorder (i.e., Anderson localization) and the strong 
correlation effects on the same footing. 

It has been known for a long time that naive mean-field approximation to effective field 
theories of localization (nonlinear <7-modcls) leads to apparent order parameter that is related 
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Fig. 1 - Typical and average DOS as a function of disorder W , for a three-dimensional cubic lattice 
at the band center (iv = 0). Results from exact numerical calculations (circles) are compared to the 
predictions of TMT (for TDOS - full line ) and CPA (for ADOS - dashed line). Finite size scaling of 
the numerical data in the critical region W = 1.17 — 1.58, and sizes L = 4 — 12 is shown in the inset, 
where pt yP (W, L)/p tyP (Wc, L) is plotted as a function of £(W)/L, and £(W) = 0.5|(W C - W)/W c \~ v 
is the correlation length in units of the lattice spacing. The numerical data are consistent with 
P = v = 1.58. 

to to the density of states, and therefore does not reflect critical behavior [H] (conceptually, 
this is due to the fact that localization manifests in transport quantities, such as the two- 
particle Green function, rather than in the one-particle Green function). This impedes, usually 
straightforward, implementation of the standard concepts of critical phenomena theory, such 
as order parameter and upper critical dimension d c . Nevertheless, the ability of the order 
parameter to label different phases, thereby delineating a complete phase diagram, makes 
it the most appealing and robust outcome of the mean-field philosophy, which survives its 
quantitative inaccuracy in dimensions below d c (there is a plethora of evidences that d c — > 
oo for Anderson localization 0). In this letter, we demonstrate how an appropriate local 
order parameter can be defined and self-consistently calculated, producing a "mean-field" 
(in the unconventional sense elaborated above) description of Anderson localization. This 
formulation is not restricted to either low temperatures or to Fermi liquid regimes. The 
promising feature of our "solution" , which we term typical medium theory (TMT) , to the 
order parameter puzzle, is in its potential usefulness outside of the realm of pure Anderson 
localization of non-interacting particles. Being amenable to easy incorporation into the well- 
known dynamical mean- field theories (DMFT) 0E] of strongly correlated systems, TMT offers 
a novel framework to address questions difficult to tackle by any alternative formulation, but 
which are of crucial importance for many electron systems of current interest dominated by 
non-perturbative physics. 

Our starting point is motivated by the original formulation of Anderson pQ , which adopts a 
local point of view |1] , and investigates the possibility for an electron to delocalize from a given 
site at large disorder. This is most easily accomplished by concentrating on the (unaveraged) 
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local density of electronic states (LDOS) pi(u>) = ^2 n S(u} — oj n )\^ n (i)\ 2 . In contrast to the 
global (averaged) density of states (ADOS) which is not critical at the Anderson transition, 
LDOS undergoes a qualitative change upon localization, as first noted in Ref. [J. This follows 
from the fact that LDOS directly measures the local amplitude of the electronic wave functions. 
When electrons localize, the local spectrum turns from a continuous to an essentially discrete 
one pQ, but the typical value of LDOS vanishes. Just on the metallic side, but very close to the 
transition, these delta-function peaks turn into long-lived resonance states and thus acquire a 
finite escape rate from a given site. According to Fermi's golden rule, this escape rate can be 
estimated pp as t~* ~ i 2 p, where t is the inter-site hopping element, and p is the density of 
local states of the immediate neighborhood of a given site. 

The typical escape rate is thus determined by the typical local density of states (TDOS), 
so that TDOS directly determines the conductivity of the electrons. This simple argument 
strongly suggests that TDOS should be recognized as an appropriate order parameter for 
Anderson localization. Because the relevant distribution function for LDOS becomes increas- 
ingly broad as the transition is approached, the desired typical value is well represented by the 
geometric average pt yp = exp{(lnp)}, where (• ■ •) represents the average over disorder. Inter- 
estingly, recent scaling analysis of the multifractal behavior of electronic wave functions 
near the Anderson transition have independently arrived at the same conclusion, identifying 
TDOS as defined by the geometric average as the viable candidate for an order parameter 
(somewhat related ideas have also been discussed earlier in Ref. L A complementary in- 
sight comes from effective field theories, applied to special models (e.g., localization on the 
Bethe lattice) that can be solved in the strong coupling limit, which point out that the full 
LDOS distribution should be considered an an "order parameter function" Despite these 
advances, a transparent technique to obtain order parameter within the framework of stan- 
dard models of localization (such as tight-binding Hamiltonian with diagonal disorder — the 
so-called Anderson model that is usually a non-interacting disordered piece of lattice Hamilto- 
nians of strongly correlated fermions) is lacking. Our principal result is shown in Fig.^ where 
comparison of the phase diagram obtained from TMT to the numerically exact treatment of 
the Anderson model on a 3D simple cubic lattice demonstrates that such approach is fully 
capable of capturing essential features (e.g., qualitative positions of the phase boundaries) of 
localization-delocalization transition occurring at strong disorder in dimensions d > 2. Thus, 
despite being an uncontrollable approximation that does not reproduce correct critical expo- 
nents in any dimensionality, TMT paves the way to include the effects of strong disorder into 
modern DMFT approaches to strongly correlated systems, where it has been a daunting task 
to go beyond coherent-potential approximation (CPA) ^3] (CPA cannot account for quantum 
interference effects induced by scattering off impurities, i.e., localization). 

To formulate a self-consistent theory for our order parameter, we follow the "cavity 
method," a general strategy that we borrow from the DMFT In this approach, a given 
site is viewed as being embedded in an effective medium characterized by a local self en- 
ergy function E(cj). For simplicity, we concentrate on a single band tight-binding model of 
nonintcracting electrons with random site energies £j with a given distribution P{si). The 
corresponding local Green function then takes the form G(u>,Ei) = [u> — Si — A(cj)]" 1 , where 
the "cavity function" is given by A(w) = A D (o; — £(cj)) with A (cj) = u> — 1/G (u>). The 
lattice Green function G (uj) = f + °° duj' D ^ } is the Hilbert transform of the bare density of 
states D(uj). Given the effective medium specified by a self-energy S(cj), we are now in the 
position to evaluate the order parameter, which we choose to be TDOS as given by 
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Fig. 2 - Typical density of states for the SC model, for disorder values W = 0, 0.25, 0.5, 0.75, 1, 1.25, 
1.275, 1.3, 1.325, 1.35. The entire band localizes for W = W c = e/2 « 1.359. 

where LDOS p(u>, £j) = — ImG(w,£j)/7r. To obey causality, the Green function corresponding 
to pty P (uj) must be specified by analytical continuation, which is performed by the Hilbert 
transform 



Finally, we close the self-consistency loop by setting the Green functions of the effective 
medium to be equal to that corresponding to the local order parameter, so that G em (uj) = 
G (u> — = Gtyp(w). It is important to emphasize that the procedure defined by these 

equations is not specific to the problem at hand. The same strategy can be used in any 
theory characterized by a local self-energy. The only requirement specific to our problem 
is the definition of TDOS as a local order parameter given by Eq. (|TJ. If we choose the 
algebraic instead of the geometric average of LDOS, our theory would reduce to CPA |14| . 
which produces excellent results for the ADOS for any value of disorder, but finds no Anderson 
transition. Thus TMT is a theory having a character very similar to CPA, with a small but 
crucial difference — the choice of the correct order parameter for Anderson localization. 

In our formulation, as in DMFT, all the information about the electronic band structure is 
contained in the choice of the bare DOS D(u>). It is not difficult to solve the above equations 
numerically, which can be efficiently done using FFT methods ■ We have done so for several 
model densities of states, and find that most of our qualitative conclusions do not depend on 
the specific choice of band structure. We illustrate these findings using a simple semicircular 
(SC) model for the bare DOS given by D{w) = - (2w) 2 , for which A D (w) = G (w)/16 [7j. 

Here and in the rest of the paper all the energies are expressed in units of the bandwidth, 
and the random site energies e\ are uniformly distributed over the interval [-W/2, W/2}. The 
evolution of TDOS as a function of W is shown in Fig. The TDOS is found to decrease and 
eventually vanish even at the band center for W = W c ~ 1.36. When W < W c , the part of 
the spectrum where TDOS remains finite corresponds to the region of extended states, and is 
found to shrink with disorder, indicating that the band tails begin to localize. The resulting 
phase diagram is presented in Fig. [21 showing the trajectories of the mobility edge (as given 
by the frequency where TDOS vanishes for a given W), and the band edge (where the ADOS 
calculated by CPA vanishes). 

Further insight in the critical behavior is obtained by noting that near W = W c it proves 
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Fig. 3 - Phase diagram for the SC model. The trajectories of the mobility edge (full line) and the 
CPA band edge (dashed line) are shown as a function the disorder strength W . 



possible to analytically solve TMT equations. Concentrating for simplicity on the band center 
(u> = 0), we can expand these equations in powers of the order parameter p a = p t yp(0) giving 
po = ap — bp 2 Q + ■ ■ •, where a — exp {—2 J deP(e) In |e|} and b = 2aP(0). The transition 
where p vanishes is found at a = 1, giving W — W c — e/2 — 1.3591, consistent with our 

numerical solution. Near the transition, to leading order p {W) — (§) 2 (W c — W), meaning 
that the order parameter exponent is = 1. The analytical solution is more difficult to 
obtain for arbitrary W. Still, the above approach can be extended to find a full frequency- 
dependent solution close to the critical value of disorder W = W c . There it assumes a simple 

scaling form pt yp (u, W) = p {W)f [uj/lo {W)], with uj {W) = J (§) (W c — W) and the scaling 



function f(x) = 1 



This is again consistent in detail with the numerical solution of Fig. [21 



(note that TDOS curves assume a simple parabolic shape close to W = W c ). 

In order to examine the quantitative accuracy of out theory, we have carried out exact 
numerical calculations for a three-dimensional cubic lattice with random site energies, using 
Green functions for an open finite sample attached to two semi-infinite clean leads ^3- We 
have computed both the average and the typical DOS at the band center as a function of 
disorder, for cubes of sizes L —4, 5, 6, 7, 8, 9, 10, 11, and 12, and averages over 1000 sample 
realizations, in order to obtain reliable data by standard finite size scaling procedures. The 
TMT and CPA equations for the same model were also solved by using the appropriate bare 
DOS (as expressed in terms of elliptic integrals), and the results are presented in Fig.Q] We 
find remarkable agreement between the numerical data and the self-consistent CPA calcula- 
tions for the ADOS, but also a surprisingly good agreement between the numerical data and 
the TMT predictions for the TDOS order parameter. For a cubic lattice, the exact value is 
W c « 16.5/12 = 1.375 /L6 , whereas TMT predicts a 20% smaller value W c w 1.1. The most 
significant discrepancies are found in the critical region, since TMT predicts the order param- 
eter exponent = 1, whereas the exact value is believed to be f» 1.58 ^7], consistent with 
our numerical data. Nevertheless, we conclude that TMT is as accurate as one can expect 
from a simple "mean- field" like formulation |19| . 

Next, we address the transport properties within TMT. The escape rate from a given site 
can be rigorously defined in terms of the cavity field A(w), and using our solution of the 
TMT equations, we find r~* = — ImA(O) ~ pt yp ~ (W c — W). To calculate the conduc- 
tivity within our local approach, we follow a strategy introduced by Girvin and Jonson |21j . 
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Fig. 4 - Conductivity as a function of the inelastic scattering rate r\ for the SC model at the band center 
and W = 0, 0.125, 0.25, 0.5, 0.75, 1, 1.25, 1.36, 1.5, 1.75, 2. The "separatrix" (a = o* independent of 
r), i. e. temperature) is found at W = W* ~ 1 (dashed line). The critical conductivity cr c (r]) ~ r] 1 ' 2 
corresponds to W = W c = 1.36 (heavy full line). 



who pointed out that close to the localization transition the conductivity can be expressed 
as a — A(yli2j42i — A11A22), where Aij = — IrnGV, is the spectral function corresponding 
to the nearest neighbor two-site cluster. We have computed a\2 by examining two sites em- 
bedded in the effective medium defined by TMT, thus allowing for localization effects. The 
vertex function A remains finite at the localization transition |21j , and thus can be computed 
within CPA |22- The resulting critical behavior of the T = conductivity follows that of 
the order parameter, a ~ pt yv ~ (W c — W), giving the conductivity exponent \x equal to the 
order parameter exponent (3, consistent with what is expected |17j . Finally, we examine the 
temperature dependence the conductivity as a function of W. Physically, the most important 
effect of finite temperatures is to introduce finite inelastic scattering due to interaction effects. 
At weak disorder, such inelastic scattering increases the resistance at higher temperatures, but 
in the localized phase it produces the opposite effect, since it suppresses interference processes 
and localization. To mimic these inelastic effects within our noninteracting calculation, we 
introduce by hand an additional scattering term in our self-energy, viz. £ — > E — ir\. The 
parameter r\ measures the inelastic scattering rate, and is generally expected to be a mono- 
tonically increasing function of temperature. The resulting behavior of the conductivity as a 
function of rj and W is presented in Fig. 21 As r\ (i.e., temperature) is reduced, we find that the 
conductivity curves "fan out" , as seen in many experiment close to the MIT . Note the 
emergence of a "separatrix" where the conductivity is temperature independent, which 
is found for W « 1, corresponding to kp£ ~ 2, consistent with some experiments 0]. At the 
MIT, wc find a c {r,) ~ p tjp (r]) ~ 77 1 / 2 . 

In summary, we have formulated a local order parameter theory for disorder-driven MITs 
that in absence of interactions reproduces most of the expected features of the Anderson tran- 
sition. In addition, the role of strong electronic correlations near disorder-driven MITs can be 
readily examined within the typical medium theory, since the local character of our approach 
offers a natural starting point for incorporating both the localization and the interaction effects 
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using the DMFT framework. 
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